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ABSTRACT 

We study the large-scale triggering of star formation in galaxies. We find 
that the largest mass-scale not stabilized by rotation, a well defined quantity in a 
rotating system and with clear dynamical meaning, strongly correlates with the 
star formation rate in a wide range of galaxies. We find that this relation can be 
explained in terms of the threshold for stability and the amount of turbulence 
allowed to sustain the system in equilibrium. Using this relation we also derived 
the observed correlation between the star formation rate and the luminosity of 
the brightest young stellar cluster 

Subject headings: galaxies: disk instabilities - galaxies: formation - star forma- 
tion: general 



1. Introduction 

Observations of normal spiral galaxies by Schmidt (1959) first suggested that their star 
formation rates (SFR) scales with the global properties. This conclusion was extended to 
other galaxies with higher SFR, such as the nuclear regions of spiral galaxies and Ultra Lu- 
minous InfraRed Galaxies (ULIRGs) by Kennicutt (1998). These galactic scale observations 
have lead to an empirical law for star formation called Kennicutt-Schmidt (KS) Law: 

Sstar OC Sgg^ , (1) 

where Sgas and Sgtar are the gas surface density and SFR per unit area. 

Since star formation is a local process that happens on subparsec scale, the correlation 
with global (galactic scale) quantities such as averaged Sgas, suggest the existence of a 
physical connection between galactic (> 1 kpc) and subparsec scales. Motivated by the 
existence of the KS law, several authors have tried to find the link in which a global/large 
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scale property of galaxies could trigger and/or regulate star formation in them (Quirk 1972; 
Wyse 1986; Larson 1987; Kennicutt 1989; Silk 1997; Tan 2000; Elmegreen 2002; Li et al. 
2005). It is also important to note that the KS law can still be explained in terms of mainly 
local processes within starforming clouds (Krumholtz et al. 2009 and references therein). 

From the study of gravitational instabilities in disks, the Toomre parameter Q = 
Csk/ttGE (or variations of that such as the star formation threshold E^rit; Kennicutt 1989) 
arises as the most natural candidate as key triggering parameter. However, the average Q 
in a galaxy is observed to be never too far from 1, from local spiral galaxies (Martin & Ken- 
nicutt 2001), to starbursts such as ULlRGs (Downcs & Solomon 1998). This is believed to 
be due to self-regulating feedback processes in the following way (Goldreich & Lyndcn Bell 
1965): if Q >> 1, then the disk will cool rapidly and form stars, while if Q << 1, then the 
star formation will be so efficient that the disk will heat up to Q ~ 1. Observed variations 
of Q (or Sgas/Scrit) are at best within factors of a few, which makes hard to explain the 
variations of 6-7 orders of magnitude that varies the SFR in galaxies, solely in terms of this 
threshold. 

However, for different disks all at the condition of marginal Toomre stabihty, is still 
possible to have different equilibrium states and this is the topic that we will explore in this 
Letter. The aim of this work, is to study why the SFR in some galaxies can be orders of 
magnitude higher than in other ones, more specifically, which galactic property triggers this 
behavior. 

This work is organized as follows. Starts with a review the gravitational instability 
analysis with a focus on the role of turbulence in §2. In §3, continues studying the effects of 
the global threshold for stability on the star formation activity, and the implications of this. 
Finally in §4, we discuss the results of this Letter. 



2. GRAVITATIONAL INSTABILITY IN REAL ISM AND THE ROLE OF 

TURBULENCE 

We start reviewing some standard results from gravitational instability analysis (Toomre 
1964; Goldreich & Lynden Bell 1965). For one of the simplest cases of a differentially rotating 
thin sheet or disk, linear stability analysis of such a system yields the dispersion relation 
for small perturbations (Binney & Tremaine 2008) of uj"^ = — 27rGS|A;| + k'^C^. Where 
Cs = \J~^ is the sound speed, S is the surface density, and k is the epicyclic frequency , 
which is given by K^{R) = i?^|l + 4^]2 (b eing Vt the angular frequency). The system becomes 
unstable when a;^ < which is equivalent to the condition Q < 1, where Q is the Toomre 
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parameter and is defined as Csk/ttGE . In such a case, there is a range of unstable length 
scales limited on small scales by thermal pressure (at the Jeans length Ajcans = C^/GT,) and 
on large scales by rotation (at the critical length set by rotation, Amt = 47r^GS/ft;^). All 
intermediate length scales are unstable, and the most rapidly growing mode has a wavelength 

2 A Jeans- 

The maximum unstable length scale in a disk, Aj-ot, is a robust quantity because it 
depends only on the surface density and epicyclic frequency of the disk and not on the 
smaller scale physics. Such scale has a characteristic mass associated, defined equals to 
Egas(Arot/2)^ and that can be expressed as: 

Mrot = r^- (2) 

AC 

Contrarily, due to the complex structure and dynamics of real interstellar medium (ISM) 
in galaxies, which cannot be described by a simple equation of state, there is not a well- 
defined Jeans length at intermediate scales. Therefore, no real lower hmit on the sizes of the 
self-gravitating structures that can form until the thermal Jeans scale is reached in molecular 
cloud cores (Escala & Larson 2008). 

One of the intermediate scales usually claimed to be a relevant in the stability of the 
ISM, is the Jeans length defined by the turbulent pressure in the medium (Aj"^^ — u^^^^^/GS). 
Besides the convenience of the idea of a turbulent pressure term that generalize the gravita- 
tional instability analysis, this term is however not a well defined quantity in ISM (Elmegreen 
& Scalo 2004). This comes from the fact that this pressure term could only be defined in the 
case where the dominant turbulent scale is much smaller than the region under consideration 
('microturbulence'), which is in fact is not the case of ISM. Rigorous analysis endeed shows 
that turbulence can be represented as a pressure only in the microturbulence case (Bonazzola 
et al. 1992). Besides turbulent motions could help to prevent collapse on some scales, those 
complex motions cannot be represented as a pressure. Therefore the gravitational instability 
analysis is not strictly applicable with a turbulent pressure term that could stabilize and 
damp all the substructure below Aj"^^. 

Although turbulence cannot be modelled as simple pressure supporting term that will 
prevent all gravitational instabilities below some scale, turbulence could prevent the growth 
of instabilities by heating the gas trough compression. However, compression from turbulence 
not only heat the gas, but also may enhance collapse and star formation on smaller scales 
(Elmegreen 2002; Krumholtz & Mckee 2005). Only when the turbulent heating becomes 
faster than the gas cooling, the net effect of turbulence is to effectively heat the gas and 
therefore prevent gravitational instabilities. Once this condition is satisfied the gas is heated 
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towards Cs ~ Vturb, the thermal pressure will stabilize and damp all the substructure below 
Aj ~ Aj"''^. If the condition Aj ~ Xrot (equivalent to Q ~ 1) is also satisfied, the system will 
then became globally stable. 

A well studied case in which both conditions could be satisfied, is in self-gravitating 
protoplanetary disks where the heating is due to gravito-driven turbulence. Gammie (2001) 

showed that turbulence heat the gas such the system becomes stable on all scales (Q > 
1), when tcooi > Stdyn, being tcooi and tdyn the cooling and orbital times. Such disk is in 
an equilibrium state (so-called ^ gravitoturbulence^) that experiences significant fluctuations, 
but the disk is stable against fragmentation and maintains itself in the brink of instability 
(Rafikov 2009) 

Galaxies are indeed observed to be close to equilibrium (Martin & Kennicutt 2001; 
Downes & Solomon 1998), with observed Toomre Q parameters never too far from 1 (av- 
eraged over the whole system and using the turbulent version of the Toomre parameter; 
Qturb = Vturb't/TrGS). This is due to self-regulation heating/cooling processes: if Qturb >> 
1, in the absence of heating driven by instabilities the disk will cool rapidly and the system 
will eventually become unstable, while if Qturb << then the self-gravity and star formation 
feedback will be so efficient that will produce enough turbulence to heat the disk towards 

Qturb ~ 1- 

However, galaxies are in an state that departs somewhat from the ' gravitoturbulenf one 
found by Gammie (2001). Because in galaxies only the turbulent Q Toomre parameter is close 
to one, not the thermal Q which is the one that guarantees stability (like in ' gravitoturbulenf 
state). Galaxies are only close to stability, the runaway growth of density fiuctuations is not 
suppressed and the formation of bound objects on different scales is ongoing (i.e. star 
formation, GMC formation, etc). They are probably oscillating around marginal stability 
due to the self-regulation feedback process, in a more dynamical fashion (and with larger 
oscillations) than the one studied in Gammie (2001). 

Besides this more complex behavior of ISM in galaxies compared with the classical 
self-regulated ^ gravitoturhulenV state studied in protoplanetary disks, the threshold is still 
well defined and there is self-regulation processes toward this marginal state. The marginal 
stability is well defined at thermal Q=l and this happens when the following two conditions 
are satisfied: Qturb ~ 1 and tcooi > theat- 
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3. Global Threshold for Stability and its Effect on the Star Formation Activity 

Besides most disks in galaxies are in the same state of being close to marginal Toomre 
stability, still some disks like nuclear disks in starbursts are able to be much more turbulent 
than others such as the disks of spiral galaxies. The reason is that is still possible to have dif- 
ferent equilibrium solutions for Q ~ 1, which is equivalent to M^ot ~ Mjeans- Simply because 
some disks starts with a higher threshold for stabihty (a larger mass-scale not stabilized by 
rotation, Mrot), their large-scale conditions requires to the self-regulation processes to pro- 
duce more turbulence in order to heat the gas towards this higher value for globally stabilize 
the system. Under this view, the self-regulation drives the system towards an equilibrium 
with a higher level of turbulence, and this more turbulent state is triggered by the initial 
condition of having a system with higher threshold for stability. 

The existence of different equilibriums in disks is particularly relevant because is believed 
that turbulence has a role in enhancing and possibly controlling star formation (Elmegreen 
2002; Krumholz & McKee 2005; Wada & Norman 2007). On its simplest form (proposed 
by Elmegreen 2002) the SFR depends on the probability distribution function (PDF) of the 
gas density produced by galactic turbulence, which appears to be lognormal in simulations 
of turbulent molecular clouds and interstellar medium (Wada & Norman 2001; Ballesteros- 
Paredes & Mac Low 2002; Padoan & Nordlund 2002; Li et al. 2003; Kravtsov 2003; Mac Low 
et al. 2005; Wada & Norman 2007; Wang & Abel 2009). The dispersion of the lognormal 
PDF, is believed to be determined by the rms Mach number of the turbulent motions if 
high-density regions are formed mainly through shock compression in a system (Vazquez- 
Semadeni 1994; Padoan et al. 1997; Nordlund & Padoan 1998; Scalo et al. 1998). Therefore 
it is expected that the SFR in galaxies scales with the velocity dispersion of turbulent 
motions. 

In summary, the existence of a higher threshold for stability (M^ot) and the self-regulation 
towards marginal stability, allows the disk to be in equilibrium with a higher turbulence 
level, which itself enhances a higher star formation activity. Therefore, both galactic-scale 
processes (Mrot) and local processes (turbulence within molecular clouds) are relevant in 
regulating star formation. This global well-defined scale in disks, Mj-ot, plays the role of 
triggering turbulence and therefore the star formation activity in galaxies. 

In order to test if this scenario is correct, we will check if the mass-scale defined by 
rotation (Eq. 1) correlates with the star formation rate as expected. For a rotationally 
supported system, the average of this mass-scale can be expressed in terms of quantities 
such as the total gas mass and gas fraction, that are easier to compare with observations 
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(Escala & Larson 2008): 

Mrot = 3 10^ Mo (^y, (3) 

where Mgas is the total gas mass in the disk and r] = Mgas/Mdyn is the ratio of the gas mass 
to the total enclosed dynamical mass within the gas radius (this varies from the disk radius 
for spiral galaxies, to the radius of the nuclear starburst disk/ring in ULIRGs). 

In figure [1^ we plot the stability threshold defined by rotation estimated from Eq. [3l 
against the measured star formation rate in those galaxies. With stars we plot the data taken 
for normal spirals (Kennicutt 1998; Kent 1987; Pisano et al. 1998; Giraud 1998; Theis 2001; 
Braine 2001; Leitherer 2002; Heifer et al. 2003; Gallagher 2005; Afanasiev 2005; Davidge 
2006; Perez- Torres & Alberdi 2007; Thilker et al. 2007), with open circles for the nuclear 
gas in normal spirals (Jogee 2005; Mauersberger et al 1996; Alonso-Herrero 2001; Hsieh et 
al 2008) and with full circles for ULIRGs (Downes & Solomon 1998). Figured^ shows a 
clear correlation between Mjot and the star formation rate, in agreement with the scenario 
outlined in this Letter. Figure [1^ supports the critical role of this threshold mass in the 
triggering of star formation. The solid line in figure [1^ represents a star formation law of 
SFR oc Mj„t 

For comparison purposes we also plot in figure [lb the measured star formation rate 
against the square of the gas fraction rj. Figure [T]d clearly shows that the scatter in this 
relation is considerably increased compared to figure [1^. This means that the gas fraction 
is not a more fundamental parameter in controlling the star formation rate. We do not plot 
the SFR against the gas mass, since it is well established that there is no correlation of the 
total gaseous mass in galaxies with their current star formation rate. 

In summary, the predicted correlation between the SFR and the maximum unstable 
mass defined by rotation is indeed observed in galaxies, in a range that spans for 5 orders 
of magnitude in SFR. This is an strong suggestion that the global threshold for instability 
indeed triggers star formation, by allowing the disk to be in equilibrium configuration with 
a higher turbulence level. 

3.1. Relation between the SFR and most luminous young stellar cluster 

Probably the most straightforward application of the scenario outlined in §2, is to check 
if its able to explain the observed correlation between the star formation rate of a galaxy 
and the luminosity of its brightest young stellar cluster (Larsen 2002; Bastian 2008). This 
is because the mass-scale studied here and that correlates with the SFR (Fig 1), also corre- 
sponds to the most massive unstable cloud in a disk, which could lead to the formation of 
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the most massive and luminous young cluster in such system (Escala & Larson 2008). 

The total luminosity of a cluster, can be computed for a given initial mass function 
(IMF) and mass- luminosity relation. Assuming a Salpeter IMF oc m~^-^^) and the usual 
mass-luminosity relation for main sequence stars (L oc m^-^), the total luminosity of a cluster 
is given by 

Ltot oc W'Z, , (4) 

where Mdoud is the total mass of the parent cloud that formed the cluster, which is assumed 
to satisfy the observed relation with the most massive star of such cluster: ~ ^cioud 

(Larson 1982). 

Taking into account the correlation found between the most massive unstable cloud and 
the star formation rate in galaxies (Fig 1), which is approximately M^^^j oc SFRt^, the 
total luminosity of the brightest cluster is given by L^^J^*^**^^* oc SFR°-^^. Finally, this can be 
expressed in terms of absolute magnitude in the V-band by My'^'^**'^* = 4.79—2.5 LogLy 
arising to the relation: 

^brightest ^_i Q^ Log SFR , (5) 

which is in good agreement with the observed slope of -1.87 (Weidner et al. 2004). Particu- 
larly good if we take into account all the uncertainties in the assumptions, such as the slopes 
of the IMF and the mass-luminosity relation. 



4. SUMMARY 

In this Letter we have studied gravitational instabilities in disks, with special attention 
on the role of turbulence in stabilizing the system. We discussed that although turbulent 
motions cannot be modelled as a pressure supporting term, turbulence can still make a 
disk globally stable when both the turbulent Toomre (Qturb >1) and Gammie (tcooi > theat) 
conditions are satisfied. This allows to the global threshold for stability, which is defined by 
largest scale in galactic disks not stabihzed by rotation, to have a clear role in the dynamics 
of the ISM. 

For a disk with a larger scale not stabilized by rotation, its dynamical equilibrium con- 
figuration (Q~l) is with a more turbulent ISM. Therefore, the role of the global threshold 
for stability in a disk is to define up to which amount the generation of turbulence is al- 
lowed. Since turbulence enhance collapse on small scales and therefore star formation, on the 
scenario proposed in this Letter, we expect a correlation between the mass-scale for global 
stability and the star formation rate. 
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We found that this relation is indeed observed in galaxies that ranges from ULIRGs to 
normal spirals. Compared with other relations such us Kennicutt-Schmidt Law, its relevance 
rely on that is the only correlation of the star formation rate with a quantity with clear 
dynamical meaning on galactic scales. 

We also explore the implications of this predicted relation between global stability and 
the star formation rate. We found that based on this relation, we are able to explain the 
observed correlation between the star formation rate in galaxies and the most luminous 
young stellar in them. 
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support from Center of Excellence in Astrophysics and Associated Technologies (PFB 06), 
from Centro de Astrofisica FONDAP 15010003 and from GEMINI-CONICYT FUND. 
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Fig. 1. — a) The star formation rate plotted against the critical mass-scale defined by rotation 
Mrot, estimated from Eq. fusing measured quantities in those galaxies. The black dots shows 
data for nuclear starburst disks, the stars for normal spirals galaxies and the open circles 
for nuclear gas in spirals. The solid line corresponds to SFR oc M.]^. b) The star formation 
rate plotted against the square of the gas fraction r] = Mgas/Mdyn, using the same symbol 
nomenclature as in a). 



